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Ij^'. 1 Introduction 

X 

With the development of observation technology, we have more opportunity to observe 
Q^. and analyze the thermal radiation from the compact star surfaces. The comparison 

Q ! between the predictions of theoretical study of the thermal evolution and observations 

of thermal emission from the compact stars gives a potentially powerful method to 
study fundamental properties of superdense matter in compact star interiors. lU |2]. 
It is generally believed that dense matter can be superfluid at such high density 
and low temperature |3]. Baryon pairing (or quark pairing) strongly suppresses not 
^ ! only the neutrino emission but also the heat capacity of nuclear matter. Meanwhile, 

f^^ I Cooper pair formation and breaking also will effect the neutron star cooling. Further- 

"^ ' more, as a neutron star spins down, some heating mechanisms may be present and 

will play important roles in the thermal evolution of NS. They have been extensively 
Q ■ discussed, for example, rotochemical heating [H |5], mutual friction between super- 

O ! fluid and normal components of the star [8], crust cracking [6] and the dissipation of 

rotational energy due to viscous damping [7] . We will focus on rotochemical heating 
in the following study. 

^ ' It's well-known ^ that for nonsuperfiuid NS with the core composed of standard 

j^ ■ nuclear matter, we have two distinct cooling regimes, slow and fast cooling. The 

transition from slow to fast cooling with increasing M occurs in a very narrow mass 
range and is sharp. After considering the superfluidity, the Urea neutrino emission 
and the NS heat capacity will be strongly suppresses, and neutron emission due to 
Cooper pairing will occur [10]. Then these effects will make NS warmer. Recently, 
when including nucleon superfluidity and using APR EOS , [TT] found a particular 
scenario of neutron star cooling that the theoretical cooling models of isolated middle- 
aged neutron stars can be divided into three distinct types: slow, moderate and fast 
cooling. They neglected the heating effects, while we will take this effects into account 
which may make the cooling scenario different. 



2 Nucleon Superfluidity 

Many microscopic studies of dense nucleon matter predict that below certain critical 
temperature nucleons will be in superfluidity, and the critical temperatures, Tcp and 
Ten, depend sensitively on the model of nucleon-nucleon interaction and many-body 
theory employed(see, e.g., [12], [IS], for references). It is widely accepted that there are 
three types nucleon superfluidities: singlet-state^ Sq pairing of neutrons (Tc = Tens) in 
the inner crust and outermost core; ^Sq proton pairing (Tc = Tcp) in the core; and 
triplet-state (^P2) neutron pairing (Tc = Tent) in the core. 

Almost all contemporary theories predict some common features of nucleon super- 
fluidity: Tcn(p) for the singlet-state neutron SF has maximum at subnuclear densities 
in the crust and vanishes at p ~ 2 x 10^^ g cm~^ while Tcn(p) for the triplet-state 
neutron SF grows up at subnuclear density, reaches maximum at p = (2 — 3) po 
(po = 2.8 X 10^^ g cm~^ is the saturation density of nuclear matter) and decreases 
with p, vanishing at p ~ (3 — 5) x 10^^ g cm~'^. Tcp also has maximum at several po 
and vanishes at higher p [lll[9|. To describe the density dependences of Ten and Tcp, 
several authors propose phenomenological models [T^ [15]. 

The nucleon superfluidity suppresses the heat capacity and neutrino processes 
involving superfluid nucleon, and what's more, superfluidity initiates an additional 
neutrino emission associated with Cooper pairing of nucleons p^ . Our codes of NS 
thermal evolution include all these effects. 



3 Rotochemical heating 

As neutron stars spin sown and contract, their structure and the chemical equilibrium 
state of beta processes change with the increasing density. Non-equilibrium reactions 
tend to restore equilibrium. During the relaxation to a new equilibrium state, the 
energy is stored. In addition, the chemical imbalance modifies the reaction rates [16j . 
If the departure from equilibrium is large enough, the net effect of the reactions is to 
increase the thermal energy at the expense of the stored chemical energy [17] . 

For npe matter in NS, their relative concentrations are adjusted by direct Urea 
reactions, 

n — )■ p + e + 17, (1) 

p + e ^ n + u, (2) 



and modified Urea reactions, 



n + N ~^p + N + e + V, (3) 

p + N + e^n + N + u. (4) 



The departure from the chemical equihbrium can be quantified by the chemical im- 
balance: 

Vnpe = SjJn " S fip - (5/ie, (5) 

where 5ni = fii—fi'^'^ is the deviation from the equilibrium chemical potential of species 
i at given pressure. In non-equilibrium state, neutrino emissivities and net reaction 
rates per unit volume of Urea processes will be modified, which can be written as [16] 



QUn,T,r^^) = Ql^{n,T)F4^), (6) 

Ar„(n,T,r^J = ^Ql\n,T)H^{^\ (7) 

where Q^'^ is the neutrino emissivity in equilibrium and rja is the chemical imbalance 
due to reaction a, T is the local temperature, n is the baryon number density, k is 
Boltzmanm's constant, and the expressions of F* and H^, are given in the appendix of 
[4j. Here we have to announce that if we consider the property of nucleon superfiuidity 
Qq'' includes the effects of nucleon superfluidity by multiplying the reduction factors 
from [9]. The total energy dissipation rate per unit volume is 

Qh = Y. ^r„r/„. (8) 



4 Thermal Evolution 

Thermal evolution of NS can be distinguished three main stages: (i) the internal re- 
laxation stage (t ~ 10-100 yr; [18j), (ii) the neutrino stage (the neutrino luminosity 
Ly ':$> L^, t ^ 10^ yr), and (iii) the photon stage {L„ <^ L^, t > 10^ yr). Af- 
ter the thermal relaxation, the redshifted temperature T°° = Te'^ becomes constant 
throughout the stellar interior. The equation of thermal evolution can be written as 

Cv'^ = -L^ -L^ + L^ (9) 

where Cy is the total stellar heat capacity, L^ is the total power released by the 
heating mechanism, L^ is the total power emitted as neutrinos, L?° is the power 
released as thermal photon. These quantities are calculated as 



L- = jdVQne^t (10) 

L~ = JdVQ^e"^, (11) 

L^ = AirR^T^'f'^ = 47^i^^(^~)^ (12) 

respectively, where dV = 4:7ir'^{gll'^)dr is the proper volume element, Q^ is the total 
neutrino emissivity contributed by reactions, a is the Stefan-Boltzman constant, R 



is the stellar coordinate radius, (ps = 0(-R), -Roo = Re'^" is the effective radius as 
measured from infinity, and T^ is the redshifted effective temperature. 

When considering the rotochemical heating effect, we can write the time evolution 
of the temperature and the chemical imbalances for NS constitute of npe matter as: 



Too = [MD{Vnpe)L^ + MM{n. 



npe 
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(13) 
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kT 



[L^HniVnpe, T) + L^jHMiVnpe, T)] + 2W^,,nn 



(14) 



where L^ = L'^^+L'^, L^^, L^, L"^ and L^ are the neutrino luminosities of neutrino 
emission processes including in our thermal evolution codes (neutrino bremsstrahlung 
in nucleon-nucleon scattering, neutrino emission due to Cooper pairing of superfluid 
nucleons, direct and modified Urea reactions) respectively, the functions M and H 
quantify the effect of reactions towards restoring chemical equilibrium, the subscripts 
D and M denote direct and modified Urea reactions, the expression of Znpe relating 
to the structure of rotating NS is in [5], the scalar Wnpe quantify the departure from 
equilibrium due to the change in the centrifugal force (oc QQ) [5l [19]. 
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Figure 1: Thermal evolution curves of 1.6Mq superfluid NS for no heating (solid line) 
or rotochemical heating with magnetic field strengths (a) B = lO^^G, (b) B = lO^^G, 
(c) B = 10"G, (d) B = 10^°G, (e) B = lO^G. The initial spin period is taken to be 
1ms. 



We plot the thermal evolution curves of I.QMq superfluid NS for no heating (solid 
line) or rotochemical heating with various magnetic field strengths (10^ — lO^'^G) in 




Log t (yr) 



Figure 2: Thermal evolution curves of I.68M0 and 1.7Mq superfluid NS with roto- 
chemical heating for different magnetic fields B = lO^'^G (dash lines), B = lO^^G 
(dot lines), B = lO^^G (dash-dot lines). The initial spin period is taken to be 1ms. 



Figure 1. And the analogs figures can be seen in [4]. However, those in [4] didn't 
include the effects of nucleon superfluidity. We can see that rotochemical heating 
increases considerably the surface temperature of stellar whose interior contains su- 
perfluid nucleon. It's quite clear that the thermal evolution curves are strongly depend 
on the magnetic field strengths, which is related with the properties of MDR. It is 
obvious that the stronger the magnetic field is the earier the interesting effects of ro- 
tochemical heating occur. For the strongest field [B = lO^^G), the spin-down is rapid 
and the heating is distinct at early times while for the lowest field {B = lO^G) the 
heating has little effect until very late times. And the most interesting heating effects 
for middle-aged stellars are those for intermediate range (10^° — lO^^G) in which the 
magnetic field of observational sources lie. 

The thermal evolution curves of I.68M0 and 1.7 Mq superfluid NSs with roto- 
chemical heating for different magnetic fields (10^^ — lO^^G) are presented in Figure 
3. We can easily see that our model is consistent with the observation data, and due 
to the different magnetic fields, the vacancy region in Figure 6 of [11] will no longer 
exist which means three distinct cooling types proposed in [11] become ambiguous. 
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Discussion 



We have studied the thermal evolution behaviors of rotating NSs with nucleon cores 
including the effects of nucleons superfluidity and the rotochemical heating. 

There are some shortcomings in our work. Firstly, we neglect processes and the heat 
capacity due to the lattice of ions in the crust. Secondly, we ignore the effect of 
superfluidity on the reaction rates. Like the neutrino emissivity, the reaction rates 
should be suppressed by superfluidity unless the departure from chemical equilibrium 
rj exceeds the sum of the energy gap of the participating baryons[20]. 



